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Abstract 



Cd i We consider a system formed by an infinite viscous liquid layer with a con- 

stant horizontal temperature gradient, and a basic nonlinear bulk velocity 



profile. In the limit of long-wavelength and large nondimensional surface ten- 
sion, we show that hydrothermal surface- wave instabilities may give rise to 
disturbances governed by the Kuramoto-Sivashinsky equation. A possible 
connection to hot-wire experiments is also discussed. 
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I. INTRODUCTION 

Thermo capillary dynamics in thin two-dimensional liquid layers, and in particular the 
convective instabilities of such flows have been a subject of much interest |]l|,3. When 
the upper free surface of a planar liquid, bounded below by a rigid plate and above by 
an interface with a passive gas, is submitted to a temperature gradient, a corresponding 
gradient in the surface tension will appear which will produce motion in the bulk fluid. If a 
vertical temperature gradient is applied to a basic static state, convective motion sets in, a 
phenomenon called Marangoni convection. However, a horizontal temperature gradient may 
also give rise to instabilities, provided the basic state is not static [Q. These are the so called 
hydrothermal instabilities, which are a coupled effect produced by both temperature and 
velocity gradients. Smith and Davis [|l],3 identified two such instabilities, which manifest 
themselves in the form of convection and surface-wave motion. 

Our concern in this letter will be the study of the hydrothermal surface-wave instability. 
It has been shown in Ref. |l[] that it is possible to have such instabilities characterized by 
a zero wave-number, provided the basic underlying flow is a nonlinear return flow. The 
long wavelenght nature of the instability allows us to broach the problem by a long-wave 
perturbative analisys, where nonlinearity may be taken into account. By using the reductive 
perturbation method of Taniuti |^, we will show that, in the limit of large nondimensional 
surface tension, the waves originated by this instability turn out to be governed by the 
Kuramoto-Sivashinsky equation @J^. This equation has already been derived in different 
physical contexts |^. In particular it has been obtained as the equation governing pertur- 
bations from a reference Poiseuille flow of a film layer on an inclined plane 0, also in the 
specific limit of large nondimensional surface tension [^. However, in this case, the insta- 
bility under consideration was hydrodynamical in nature, and not hydrothermal, as is the 
case considered in the present work. 

The basic interest of the results presented in this work is concerned with recent ex- 
perimental evidences for the existence of surface hydrothermal waves. For example, it has 



been reported in the literature the observation of such waves in a cyhndrical container . 
Moreover, hot-wire experiments performed recently pO| , p!l| have indicated the presence of 
propagative patterns. In the final section, we will speculate on a possible connection between 
our results and these experiments. 

II. MATHEMATICAL FORMULATION 

We consider a fluid layer of height d, bounded below, at 2; = 0, by a rigid perfectly 
insulating plate, and above, ai y = d + r]{x, t), by a free deformable surface in contact with 
a passive gas of negligible density and viscosity. The liquid is characterized by a density 
p, thermal conductivity k, thermal diffusivity k, unit thermal surface conductance h, and 
dynamic viscosity p. To the free surface we associate a surface tension T, which will be 
assumed to depend on the local temperature 9 according to the linear law 

T = To-^{e-9o), (1) 

where 7 is a positive constant, and Tq, 60 are reference values for surface tension and tem- 
perature, respectively. 

We will be concerned with effects coming from thermocapillarity only. Therefore we 
will neglect gravity. This is a good approximation for a thin enough layer, or a layer in a 
microgravity environment. The equations governing the fluid motion are written as: 

Ux + Wz = 0, (2) 

p {ut + UU.J, + WU:,) = -p^ + p {u^x + u^z) , (3) 

p (Wt + UW:^ + WWz) = -Pz + /i (W^x + Wzz) , (4) 

9t + u9^ + w9z = K {9^^ + 9^^) , (5) 



where the subscripts denote partial differentiation, u and w are, respectively, the x and z 
components of the velocity, and p is the pressure. On the upper free surface, these equations 
are subject to the following boundary conditions 
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Equation (6) is the kinematic boundary condition, Eqs. (7) and (8) represent the normal 
and tangential stress balance at the upper free surface, and Eq. (9) is the equality of the 
heat flux at the interface, where 9oo is the temperature aX z -^ oo. On the lower plate, the 
boundary conditions are 



u = w = 6^ = Q 



(10) 



which characterizes the no slip and zero heat flux conditions. 

Perturbation theory will be performed on a basic state of the system, defined by imposing 
to the infinite liquid layer a constant horizontal temperature gradient 



dx 



(11) 



with (3 a constant, which will give rise to a basic velocity profile in the bulk fluid. Moreover, 
we will assume a nonlinear velocity profile since, as remarked in Ref . |1| , a linear flow profile 
is not susceptible to long-wave instabilities. 

We now proceed to a nondimensionalization of the equations and boundary conditions. 
To this end, we take d as the unit of lenght, /i/7/3 as the unit of time, (3d as the unit of 



temperature, 7/3 as the unit of pressure, and Tq as the unit of surface tension. In this process 

the following dimensionless numbers appear: 

P7/3rf2 II hd pdTo 

R = — , cr = — , B = —- , S 



/X^ pK k jl'^ 

where R is the Reynolds number, a is the Prandtl number, B is the Biot number, and S is 
the surface tension number. In addition, we have also the Marangoni number M, which is 
defined hy M = Ra. 

The basic state, from which we will consider small disturbances, is the return flow solution 
to the Eqs.(2-10) which, in the limit S —>■ oo, is given by: 

(12a) 
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Ooo = -x] Vb = 0. (12d) 

This solution has been obtained in Ref. |12|. It is exact in the limit of infinite surface- 



tension. However, in what follows we will consider 5* to be large, albeit finite. Accordingly, 
we should introduce corrections to Eqs.(12) up to some appropriate order of S^^. These 
corrections, though, are of higher orders than those contributing to our calculations, except 
in Eq.(7). In this equation, the right hand side is proportional to S. To match orders 



correctly we make use of the basic solution i^b correct up to order S , which is given by fT^ 



R n R 

^^ = -4S" +165"- 
It is important, at this point, to note that the range of validity of the basic solution is not 
arbitrary. As shown in Ref. [Q , the largeness of 5* is connected to the largeness of the aspect 
ratio A, the ratio of the width of the layer to its depth. In fact, we must have S ~ 0{A'^) 
as A — > oo, so that we are not allowed to take x ^ oo independently of S. 



III. PERTURBATION THEORY 

Let us consider the following perturbations to the basic state : 

u = Ub + e{uo + eui + ....) , (13a) 

w = e^ {wo + ewi + ...) , (13b) 

P = Pfe + e(po + epi + •••) ) (13c) 

e = 9h + e{9o + e9i + ...) , (13d) 

r] = r]b + e{r]o + er]i + ...) , (13e) 

where e is a small parameter. Next, we introduce the slow variables according to 

^ = e{x-ct), (14) 

T = eH, (15) 

and we suppose that the perturbations depend on (x, t) through (^, r) only. We now make 
more precise the meaning of large albeit finite for S, by assuming that 

S = e'^S, (16) 

with S r^O{e^). 

We are able, at this point, to obtain an order by order solution to the problem by 
integrations in the z variable. At each order, the bulk equations and the boundary conditions 
at the bottom yield u, w, 6,p in terms of three arbitrary functions, which will be determined 
in terms of rjo by the boundary conditions at the upper free surface. However, as we have 
four equations, a compatibility condition will arise at each order. Explicitly, at order e°, we 
get: 
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From Eqs. (6) and (8), we get as the compatibility condition at this order 
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At the next order, the relevant solutions are 
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with f = f{C,, t) an arbitrary function. Also at this order, from Eq.(6), we have : 
3 1, 5 R S 



From Eq.(8), we get: 
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4^iC + 2/^ = j(^+2J^o,, + ^r,o,«,. 



Substituting into Eq.(^), we obtain an evolution equation for rjo: 
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Returning to the laboratory coordinates (x, t), and remembering that S = e^S and rjo = e ^rj, 
Eq.(|2D reads 



Vt+ [C- -VJVx + R [j^ + J^j Vxx + ^Vxxxx = , 



(23) 



Defining a new field by 
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By a simple rescaling of variables, this equation can be written as 



(24) 



Ct + CCx + Cxx + Cx 



0. 



(25) 



which is the Kuramoto-Sivashinsky (KS) equation in its standard form 

The linear analysis performed by Smith and Davis [Q can be recovered from Eq.(24) if 
we neglect the nonlinear term, and take 



( ~ exp [i {ax — cut)] . 



This leads to the dispersion relation 
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from which we can see that the second derivative and the fourth derivative terms of Eq.(p4D 
are, respectively, antidissipative and dissipative. Therefore there exists a critical point at 
the long wavelenght limit, given by cj = 0, in which these terms compensate each other. In 
this limit, a — ;■ 0, and assuming that S = a~^S, we obtain the critical Reynolds number 
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which is the result of Smith and Davis ilH . 



IV. DISCUSSION 



The KS equation is a non-integrable equation, despite the fact that some exact solutions 
have been found PJTB[|. Generically, the solutions to the KS equation are characterized 
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by the coexistence of coherent spatial structures with temporal chaos |T^. In particular, 
the search for a better understanding of the long wavelength properties of the KS equation 
has been intensively pursued lately [|1^]. And it is precisely in the long- wave limit that we 
have obtained the KS equation as governing hydrothermal surface-wave instabilities in a 
thin viscous layer of fluid with horizontal temperature and velocity gradients, and in the 
limit of large nondimensional surface-tension. This result could possibly be connected to a 
recent experiment |^ which has been performed to study the features of hydrothermal waves 
appearing when a cylindrical liquid layer is laterally heated. The existence of hydrothermal 
wave instabilities has been observed experimentally, and the results were found to partially 
agree with the previous theoretical predictions of Smith and Davis |jl[. If we believe that, 
despite the difference in the geometry the results presented in Ref. can somehow be 
considered as an approximation to the theoretical predictions, then the long-wave dynamics 
of the free surface described by the KS equation should play a role in the explanation of the 
propagative patterns observed in this kind of experiment. 

At last, we would like to speculate on a possible connection between the results obtained 
in this letter and hot-wire experiments |]TD|,|TT|. These experiments consist of a hot wire 
placed horizontally just below a free surface of a fluid. Above a certain critical value of the 
electrical power supplied to the hot wire, propagating patterns are observed. The existence 
of both horizontal temperature gradients and a basic velocity flow, as the fluid is convecting, 
makes contact with the theoretical setting in which hydrothermal waves show up. Of course, 
this is not enough to ensure that the above results can have an immediate application to 
the hot-wire experiments. However, since the system under consideration here presents 
similarities with the above mentioned experiments, they could help to elucidate the physical 
mechanism responsible for the observed propagative patterns. The idea to mimic features of 
a convecting fluid by a two-dimensional system with a velocity gradient has been examined 
by Roze |T^, based on data of an experiment with a much shorter wire |T^. However, this 



approach, which used a linear flow profile and linearized equations, turned out to be rather 
disappointing. In fact, as shown in Ref. |jl|, there is no long-wave instability in this case. 

9 



Only when a nonlinear return flow solution is supposed, which is indeed more appropriate 
for a closed system, long-wave hydrothermal phenomena appear. Therefore, we think it 



would be interesting to reexamine the ideas proposed in Ref. [0 for the hot wire problem 
in the light of the present results. 

In summary, we have made in this letter a nonlinear analysis of the long wavelength hy- 
drothermal instabilities, which can be considered as an extension of a problem flrst discussed 
by Smith and Davis [jl| in a linear approximation. We have shown that the free-surface dy- 
namics is governed by the Kuramoto-Sivashinsky equation, and that the previous results 
can be recovered from ours as a particular case. 
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